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Abstract— The  step  response  characteristics  of  first  and 
second  order  systems  are  well  known.  On  the  other  hand , 
the  step  response  of  fractional  order  systems  (FOSs)  with  2- 
termfractionaldenominator  is  like  those  of  first  and  second 
order  systems.  But  there  are  important  differences  between 
the  two  types  of  characteristics.  Considering  the  step 
response,  the  behavior  of  simple  FOS  with  a  denominator 
polynomial  having  unity  term  and  the  other  involves 
fractional  power  is  investigated  in  this  paper  comparatively 
with  1st  and  2nd  order  systems.  The  results  bring  light  for 
the  design  of  fractional  order  control  systems  (FOCSs). 
Keywords —  Fractional  order ,  control  system ,  step 
response,  time  constant,  rise  time,  overshoot. 

I.  INTRODUCTION 

With  increasing  technology  and  high-speed  computers,  new 
numerical  methods  for  modelling  and  solving  physical 
systems  have  been  presented.  FOSs  have  founded  many 
applications  not  only  in  this  respect,  but  they  also  constitute 
more  realistic  models  for  physical  phoneme.  Therefore, 
plenty  of  literature  has  appeared  for  analyzing  and 
designing  these  systems  in  the  last  two  decades  [1-5]. 
Especially,  fractional  order  proportional  integral  derivative 
(FOPID)  controllers  have  appeared  extensively  in  control 
system  design  [6-8]. 

As  some  examples  on  very  recent  literature,  [9]  presents  a 
FO  positive  position  feedback  compensator  and  it  is 
compared  to  the  standard  integer-order  one  for  active 
vibration  Control  of  a  rectangular  free-edged  carbon  fiber 
composite  plate.  In  [10],  FO  control  of  dissimilar  redundant 
actuating  system  is  used  in  large  passenger  aircraft;  two 
dissimilar  actuators  used  for  position  control  surface  are 
modelled  by  fractional  orders  and  three  FO  controllers  are 
used  to  improve  transient  response  of  the  system.  In  [11],  it 
is  aimed  to  solve  the  stability  problem  for  some  FO 
nonautonomous  systems  by  linear  state  feedback  control 
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and  adaptive  control;  a  new  property  for  Caputo  fractional 
derivative  is  used  to  derive  some  sufficient  conditions  for 
the  global  asymptotical  stabilization. 

Considering  sophisticated  and  rather  special  applications 
[12-17],  there  are  a  huge  number  of  publications  about 
FOSs ;  so  only  very  few  of  them  are  referred  in  this  paper.  On 
the  other  hand,  a  compact  publication  investigating  the  step 
response  characteristics  of  even  simple  FOSs,  namely  rise 
time,  settling  time,  delay  time,  overshoot,  oscillation  period, 
damping  time  constant  of  oscillations,  and  some  others 
which  do  not  exhibit  similarity  with  those  of  integer  order 
systems  [18]  is  not  yet  present.  The  purpose  of  this  paper  is 
to  draw  attention  to  the  differences  between  the  step 
response  characteristics  of  simple  FOSs  and  1st  and  2nd 
order  systems  and  thus  to  supply  some  lights  for 
FOCSs’  designers. 

For  the  mentioned  purpose,  time  domain  characteristics  of 
first  and  second  order  systems  are  reviewed  in  Section  2.  In 
Section  3,  the  step  response  characteristics  of  o  a  FOS 
with  2-term  denominator  are  presented  depending  on  the 
fractional  power.  Section  4  discusses  and  lists  the 
differences  between  2nd  order  and  fractional  order 
characteristics.  Section  5  covers  the  conclusions. 


II.  REVIEW  OF  FIRST  AND  SECOND  INTEGER 
ORDER  SYSTEMS 

It  is  well  known  that  the  simple  first  order  transfer  function 

h'(s)=At  (la) 

has  a  step  response 

y1  (t)  =  1  -  e~t/Vl  (lb) 

monotonically  increasing  from  0  to  1  with  a  time  constant 
t  =  Pi  [18].  The  dc  gain  is  assumed  to  be  1  which  is  not 
important  for  relative  time  domain  characteristics. 

For  the  second  order  transfer  function, 
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H2(s)  =■ 


9  ^  — ;  (2a) 

5  +  2f(DnS  +  0)^ 

is  the  standard  form  where  o)n  is  undamped  natural 

frequency  and  f  >  0  is  the  damping  ratio  [18]. 

Overdamped  step  response 


1  /  __L  __L\ 

y(t)  =  1  — - (r2e  T2  -  TlJ,  (2b) 

T2  -  Ti  V  / 

occurs  for  for  f  >1,  where  two  time  constants  are  = 
+  ^2-l)andr2  =  l/<wn(f  -  Vf2  -  l)  so  that 
it  increases  monotonically  from  0  to  its  final  value  1.  For  £ 
=  1,  the  response  is  similar  in  shape  to  the  overdamped 
case,  but  the  variation  is  described  by 

y(t)  =  1  -  e_"nt(o)n  +  l)  (2c) 

which  is  known  as  critically  damped  case. 

Underdamped  response  occurs  for  0  <  ^  <  1.  The  response 


p  S  'JJ7TL  s  r - 

y(t)  =  1  -j^Tsin(y/l-$2cont  + 


sin  1VT^)(2d) 
is  stable  and  approaches  to  the  reference  value  1  with 
damped  oscillations  with  frequency  yjl  —  <f2a)n  and 
damping  time  constant  l/fo)n.  The  special  case  of 
oscillatory  response  occurs  for  =  0,  the  system  is 
undamped,  and  the  step  response  is 

y(t)  =  1  -  coscont  (2e) 

which  represents  sustained  oscillations  with  frequency  con 
and  undamped  oscillation  period  Tn  =  2n  /  con. 

All  the  step  responses  considered  so  far  stars  from 0  at  time 
t  =  0  and  approaches  to  1  monotonically  or  oscillatory 
manner  as  Lim  t  ->  oo  (except  the  undamped  case  =  0). 
The  following  time  domain  characteristics  are  defined  for 
these  responses: 

Time  constant  r:  Time  required  for  the  response  to 
reach  1  —  ( 1/e )  of  its  final  value. 

Rise  time  Tr:  Time  required  for  the  response  to  reach 
from  10  %  to  90  %  of  its  final  value. 

Settling  time  Ts:  Time  required  for  the  response  to 
stay  around  its  final  value  with  an  error  less  than  2  %. 

Oscillation  period  T0\  Period  of  oscillations  for 
underdamped  case. 

Peak  time  Tmax\  The  time  the  response  reaches  its 
maximum  value. 

Peak  value  vmcix:  Maximum  value  of  the  response; 

y(Tmax)  • 

Overshoot  vcrh  :  ymax  —  1 ;  How  much  the  response 
exceeds  the  final  reference  1. 

Percent  overshoot  POSH:  It  is  defined  by  {yosh/ 
Vref) 100  • 
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Reduction  Ratio:  RR:  Ratio  of  successive  overshoots . 
From  Eq.  (lb),  it  is  straightforward  to  show  that 


II 

(3a) 

Tr  =  ln9p±  =  2.197225  plf 

(3b) 

Ts  =  ln50p1  =  3.912023  p±. 

(3c) 

In  summary  for  a  first  order  system;  i)  The  step  response 
increases  exponentially  to  its  steady  state  value  without 
any  oscillations,  ii)  Rise  time  and  settling  time  are  some 
multiples  of  time  constant  t  =  p1. 

From  Eq.  (2b),  the  over  damped  response  of  a  second  order 
systems  has 

r :e~T(°n^ sinh  [ra)nA/<f2  -  1  +  cos/i-1<f] 

=  e-y<f2-l.  (4a) 

Tr  —  t2  ~  tltwhere  (4b) 

sinh[tl(Jon^2  -  1  +  cos/i-1<f]  =  0.9 -  1, 
e-t20)^sinh\t2(onyf^  -  1  +  cosh -1f]  =  0.1^/f2  -  1 
Ts\  e~TsC°n^ sink  \rscon^2  —  1  +  cos/i-1f] 


=  0.02/f2  -  1.  (4c) 

From  Eq.  (2c),  the  critically  damped  response  of  a  second 
order  systems  has 

t  =  [l  +  ln(l  +  a)n)]/mn,  (5a) 

Tr  =  ln9  /con,Ts  =  ln[50(l  +  a)n)]/a)n,  (5b,  c) 


T  = 

1  s 


-In  ■ 


50 


■  ln- 


50 


fan”  J  1-f2  (5d) 

From  Eq.  (2d),  the  under  damped  response  of  a  second 
order  systems  has 

t:  e~TC°n^ sin  [ro)nA/l  -  <f2  +  sin'll  -  <f2] 

=  e-yi-f2, 

Tr  —t2  —  tlf  where 

e-tl(°nZsin  [t16c)riA/ 1  —  ^2  +  sin'll  -  f2] 

=  0.9^1  —  (2,e"t2&,^sm  [t2^nV^-  “ 

+  5m_yi  —  <f2]  =  O.l^/l  —  <f2, 

50  2n 

T 

>  1  o 


(6a) 

(6b) 


7)  = 


-In 


fan  Vl  -f2  Vl 


T  =  ■ 

1max 


CO- 


’nVW^ 


>  Vrnax  1  "b 


0)„ 

ziU 

Vi-f2 


Vi-f2 


(6c,  d) 

(6e,f) 


gVT^F  lOOe^-f2 

a/W7 

-2f7T 

/?/?  =  (6i) 

For  the  succeeding  peaks,  Eqs.  (6e-i),  should  be  modified 
by  using  (2  k  —  l)n  instead  of  n  where  k  represents  the 
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peak  numbers;  for  the  first  peak  k  =  1  and  Eqs.  (6e-i)  are 
valid. 

From  Eq.  (2e),  the  un-damped  response  of  a  second  order 
systems  has 


T 

ir 


(x)n 

cos-1(0.l)  —  cos-1(0.9) 


(7a) 

(7b) 


T  =  oo  T  —  —  T  —  — 

1  s  >  1  o  .>  1  max  . ,  > 

a)n  0)n 


(7c;  d,  e) 


ymax  yosh  1,  (7f,  g) 

POSH  =  100,  RR  =  1.  (7h,i) 

To  be  able  to  make  comparisons  with  the  step  response 


characteristics  of  the  fractional  order  system  dealt  with  in 
the  next  section,  and  as  well  as  for  the  sake  of 
completeness,  the  step  responses  of  the  second  order 
system  in  Eq.  (2a)  are  shown  in  Fig.  1  for  a  few  values  of 
<f.  In  the  same  respect,  some  of  the  time  domain 
characteristics  are  plotted  in  Fig.  2.  The  derived 
transcendental  equations  for  t  and  Tr  (see  Eqs.  (4a, 6a)  and 
Eqs.  (4a, 6b),  respectively)  are  solved  numerically  to  obtain 
the  plots  for  relevant  characteristics  in  Fig.  2.  This 
approach  is  not  needed  for  obtaining  the  remaining 
characteristics  in  the  figure  for  the  necessary  formulas  are 
presented  explicitly. 


Fig.l :  Step  responses  of  the  system  in  Eq.  (2a)  for  values  off;  =2, 1,0.2,  0;  a)n  =  1. 


Fig.2:  Some  time  domain  characteristics  of  the  systemin  Eq.  (2a)  for  f  G  [0,3]. 
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III.  INVESTIGATION  OF  THE  SIMPLE 
FRACTIONAL  ORDER  SYSTEM 

Consider  the  fractional  order  transfer  function  with  a 
constant  numerator  and  2-term  fractional  denominator 

=  prr?  («) 

When  a  =  1  (2),  the  transfer  function  is  a  1st  order  (2nd 
order  with  con  =  1,  =  0)  integer  type  system;  hence, 

the  step  response  is  an  increasing  exponential  (sustained 
oscillation).  But  there  are  differences  in  the  step 
responses  for  the  general  values  of  a: 

Fig.  3  shows  the  step  responses  for  a  E  [0.01,1.99].  All 
the  responses  tend  to  the  reference  value  1  as  t  ->  oo, 
since  inserting  5  =  0  and  assuming  a  =f=  0  in  Eq.  (8),  dc 
gain  of  this  system  is  seen  as  1.  For  small  values  of  a 
(for  eg.  a  =0.01),  the  step  response  is  almost  equal  to  V2 
(the  transfer  function  is  equal  to  a  constant  gain  of  V2  for 
a  =  0);  except  t  =  0  and  t  ->  00;  for  these  values,  the 


response  starts  from  0  and  tends  to  the  reference  value  1 
(this  can  not  be  seen  in  the  figure  since  the  final  value  of 
t  is  chosen  30  in  simulations).  In  the  figure,  exponential 
like  increase  of  step  response  for  0  <  a  <  1,  exponential¬ 
like  decaying  oscillations  for  large  a  (for  eq.  a  =  1.99) 
are  obviously  seen.  For  a  >  1  but  not  too  much,  there  is 
an  overshoot  in  the  response  then  after  exponential -like 
decaying  to  1  without  oscillations  are  also  seen.  Fora  > 
2,  the  system  is  unstable,  and  the  step  response  increases 
exponentially  and  oscillatory  like  manner,  which  is  a  case 
not  shown  in  the  figure. 

We  note  that  those  responses  for  1  <  a  <  2  and  0  <  a  < 
1,  although  resemble  to  those  of  a  second  order  and  of  a 
first  order  integer  order  systems,  respectively,  there  are 
important  differences. 


1.8 


1.6 


1 A 


,  1-2 


0.6 


15 

Time 


Fig. 3:  Step  responses  of  the  FOS  in  Eq.  (8)  for  different  values  of  a. 


For  example,  if  there  exist  oscillations,  the  oscillation 
period  is  not  constant  and  naturally  the  peak  times  do  not 
occur  uniformly  as  shown  in  Fig.  4.  This  figure  shows  the 
change  of  oscillation  periods  starting  from  each  Tmax. 
Numerical  data  show  that  the  first  oscillatory  like 
response  (number  of  maxi  mu  ms  =  2)  occurs  for  a  — 
1.34,  for  which  the  period  of  oscillation  is  6.1750.  But 
this  case  is  not  observed  on  the  graph.  The  plot  starts 
from  the  case  where  there  exists  3  maximums  and 


resultantly  2  oscillation  periods;  this  case  occurs  fora  = 
1.49  for  which  the  first  period  is  7.188  and  the  second 
one  is  5.4470.  The  plot  for  this  case  is  indicated  in  the 
figure.  It  is  seen  that  the  first  oscillation  period  for  each 
a  is  maximum  and  it  decreases  with  increasing  a.  Further, 
for  each  a  the  oscillation  period  decreases  for  subsequent 
oscillations.  This  decrease  reaches  to  its  highest  value  just 
before  the  oscillations  stop.  For  a  =  1.99  =  2,  period  of 
oscillations  is  almost  constant  being  6.284  for  first 
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oscillation,  6.283,6.283,6.284  for  the  following  3  though  expected  since  a  =  1.99,  may  originate  from 

oscillations.  The  change  in  the  3rd  digit  after  decimal,  numerical  errors. 


Fig. 4:  Change  of  oscillation  period  with  Tmax. 


The  exponential  decays  or  increases  do  not  have  fixed 
time  constants;  hence  reduction  ratio  decreases  for 
succeeding  maximum  values  as  seen  in  Fig.  5.  This  figure 
shows  the  reduction  ratio  of  subsequent  peaks.  Reduction 
ratio  starts  from  a  =  1.34  where  at  least  2  peaks 
(maximums)  occur  for  the  first  time  (Reduction  ratio  = 
14.4192)  and  ends  at  a  =  1.99  where  5  peaks  occur 
until  t  =  30  s.  Hence,  reduction  ratio  is  defined  for  the 
first  4  peaks.  In  the  figure,  the  plot  for  a  = 


1.34, 1.35,  ...,1.48,  reduction  ratio  can’t  be  plotted  (since 
point  plot  is  not  used).  Therefore,  the  first  curve  starts 
from  a  =  1.49  for  which  3  peaks  occur,  and  2  reduction 
ratios  are  defined  (20.0430  and  3.3023).  It  is  seen  that 
the  reduction  ratio  from  the  first  peak  to  the  second  one 
decreases  as  a  is  changing  from  1.49  to  1.99.  Further,  the 
ratio  of  the  subsequent  peaks  for  each  a  decreases  with 
time  (or  the  oscillation  number). 
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The  size  of  overshoots  is  shown  in  Fig.  6.  For  a  =  1.2, 
there  is  a  single  overshoot  which  is  0.0744;  this  is  not 
shown  in  the  graph  since  there  is  not  a  second  one  so  that 
response  is  not  oscillatory. 

In  Fig.  6,  the  values  of  overshoots  are  plotted  against  the 
time  of  maxi  mu  ms.  The  first  maximum  occurs  for  a  — 
1.01 ;  since  a  single  point  plot  is  not  shown,  the  plot  in  the 
figure  is  started  from  a  =  1.34  where  the  second  peak 
occurs  for  the  first  time;  (for  a  =  1.34,  the  first  overshoot 
and  the  second  one  are  0.1640  and  0.0114,  respectively). 
Hence,  for  a  G  [1.01,1.33]  there  is  one  maximum  in  the 
step  response  and  no  other  maximums  occur  until  a  = 
1.34.  Note  that  this  last  statement  is  confined  to  the 


present  numerical  data.  More  elaborate  numerical 
simulation  covering  values  of  a  to  four  decimal  digits 
show  that  this  interval  is  a  G  [1.0001,1.3395],  and  no 
other  maximums  occur  until  a  =  1.3396. 

In  the  simulation  results  shown  in  Fig.  6,  for  a  G 
[1.63,1.99],  there  are  five  peaks  in  the  step  response.  The 
overshoots,  as  a  general  rule,  change  so  that  for  each  a 
subsequent  peaks  decrease,  while  the  rate  of  decrease  gets 
stronger.  Note  that  among  the  peaks  occurring  for  a  = 
1.67,  the  minimum  one  is  the  fifth  one  and  it  is  equal  to 
0.0001018660;  this  point  could  not  be  data  tipped  due  to 
the  crowdedness  of  curves  near  to  it,  the  minimum  peak 
value  is  found  from  the  workspace  data. 


Fig.  7  shows  the  maximum  value  of  step  response  against 
its  occurrence  time.  This  figure  is  plotted  for  the  case  of 
responses  that  have  at  least- 2- maximum  values  for  t  G 
[0,  30].  These  responses  are  for  a  =  1.34,1.35,  ...,1.99. 
a  =  1.34  is  the  first  value  of  a  where  at  least  2 
maximums  (hence  oscillatory  like  response)  occur.  More 
elaborate  numerical  simulation  shows  that  this  value  of  a 
is  a  =  1.3396  (as  mentioned)  for  which  oscillatory 


damping  starts  with  period  of  oscillation  6.0680.  After 
the  oscillations  start  at  a  =  1.34,  the  number  of 
oscillations  increases  with  increasing  a  so  that  5 
maximums  occur  for  a  =  1.99  in  [0, 30].  The  time  of 
occurrence  of  the  maximums  are  not  the  same;  they  get 
near  to  each  other  as  a  ->  1.99,  and  they  get  more 
different  values  as  a  ->  1.34. 
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Fig. 7:  Maximum  value  of  step  response  against  the  time  of  occurrence. 


Next,  it  is  worth  to  observe  the  change  of  first  oscillation 
period  and  the  number  of  peak  values  in  the  first  30  s. 
These  are  plotted  in  Fig.  8 .  It  is  observed  that  there  are 
no  oscillations  in  the  step  response  for  a  =  0,  0.5, 1.  For 
a  =  1.2,  there  are  no  oscillations  but  there  exist  1 
maximum  with  an  overshoot  of  0.0744  as  mentioned 


before.  For  a  =  1.7,  a  =  2,  and  a  =  2.018  there  are  5 
maximums  in  the  interval  tG  [0,30];  whilst  period  of 
first  oscillation  changes  with  a,  i.e.,  it  is  6.543  for  a  = 
1.7,  and  it  is  6.283  for  both  a  =  2  and  a  =  2.018.  Note 
also  that  period  of  succeeding  oscillations  is  not  constant 
for  the  same  a. 


Fig.  8:  Period  of  first  oscillation  and  number  of  peaks  in  30  s  against  a. 

In  Fig.  9,  the  variation  of  the  time  constant  r,  rise  time  Tr,  and  settling  time  Ts  are  plotted  for  a  =  0.5, 1, 1.2, 1.7.  It 
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is  seen  that  r  is  the  least  effected  by  a.  Tr  could  not  be  time  gets  large  very  rapidly  with  increasing  a\  infact  for 

recorded  for  a  =  0.5  since  it  is  greater  than  total  run  time  a  >  2  it  is  oo.  Therefore,  the  data  for  a  =  2  and  a  = 

t  =  30  s.  For  a  =  1,  r  <  Tr  <  Ts,  with  values  1.001  <  2.018  is  not  included  in  the  plot. 

2.198  <  3.913.  It  is  important  to  note  that  the  settling 


Fig.9:  Variation  of  time  constant  (t  =  Tau  ),  risetime  (Tr)  and  settling  time  (Ts)  with  a  included  in  the  plot. 


Finally,  in  Fig.  10,  the  damping  ratio  of  a  2nd  order  figure  can  be  used  for  finding  the  values  of  f  and  a  for 

system  and  the  fractional  power  a  of  the  fractional  order  the  desired  overshoot, 

system  are  compared  with  respect  to  the  overshoot.  The 


The  step  response  characteristic  investigated  so  far  are 
already  enough  to  highlight  the  important  properties  and 
differences  between  some  simple  integer  order  and 
fractional  order  systems.  Therefore,  this  section  is  ended 
at  this  step.  However,  for  a  continuous  variation  of  step 


response  characteristics,  such  as  period  of  oscillations, 
number  of  peaks,  t,  Tr,  Ts,  POSH  with  the  system 
parameters  a  and  £,  as  well  as  variation  of  a  for  fractional 
order,  and  f  for  second  order  systems  with  the  overshoot  can 
be  obtained,  these  are  not  included  for  the  length 
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limitations  of  the  paper.  Further,  they  are  not  vital  to 
observe  the  following  comparisons. 

IV.  DIFFERENCES  BETWEEN  2ND  ORDER  AND 
SIMPLE  FRACTIONAL  ORDER 
CHARACTERISTICS 

We  have  the  apparent  observation  from  the  previous  two 
sections  that  the  following  important  similarities  and 
differences  appear  between  the  step  responses  of  a  second 
order  system  and  the  simple  FOS  considered.  The 
comparison  is  done  on  the  base  of  the  transfer  functions 
in  Eq.  (2a)  for  a  second  order  integer  system  and  Eq.  (8) 
for  a  FOS.  Being  considered  before,  the  extreme  cases 
a  =  0,  a  =  0  ,  and  a  =  2  are  not  included  in  the 
comparison. 

1)  All  the  time  domain  characteristics  of  a  second  order 
system  are  express ed  by  analytic  formulas  as  given 
in  Eqs.  (4-7).  For  FOS,  there  does  not  exist  such 
analytical  formulas;  only  graphical  dependences  on 
the  fractional  power  a  are  available.  See  Figs.  4-9. 

2)  For  1  <  and  0  <  a  <  1,  both  of  the  systems 
exhibit  monotonically  increasing  step  responses 
starting  from  0  and  tending  to  the  final  value  1  as 
Lim  t  -»  oo.  There  is  no  overshoot  in  the  responses. 
See  Fig.  1  for  2nd  order  system,  and  Fig.  2  for 
fractional  one. 

3)  ForO  <  <  1  and  1  <  a  <  2,  both  systems  generate 

overshoots  which  decay  to  zero  andthe  step  responses 
tend  to  the  final  value  1  as  Lim  t  ^  oo.  See  Fig.  1 
and  Fig.  2,  respectively. 

4)  For  0  <  £  <  1  and  1  <  a  <  2,  although  the  step 
responses  of  both  the  second  order  system  and  that 
of  the  FOS  are  similar  in  the  sense  that  they  are  not 
monotonically  increasing,  they  have  overs  hoots,  and 
they  approach  to  unity  at  infinite  time,  they  have 
completely  different  characteristics: 

a)  There  are  infinitely  many  overshoots  in  the 
step  response  of  the  2nd  order  system,  whilst 
the  number  of  overshoots  in  the  step  response 
of  the  fractional  system  increases  with  a.  For 
example,  there  is  only  one  overshoot  and 
naturally  one  maximum  for  1  <  a  <  1.3395, 
and  the  second  overshoot  starts  for  a  = 
1.3396 .  Refer  to  Fig.  8  to  see  how  the  number 
ofmaximums  (oscillations)  increases  with  a. 

b)  The  period  of  sustained  oscillations  and  their 

damping  rate,  that  is  RR,  are  constant  for  each 
0  <  <  1  and  they  do  not  change  with  time; 


they  are  given  by  Eqs.  (6d)  and  (6i), 
respectively.  On  the  other  hand,  these 
quantities  change  with  time  for  FOS;  see  Fig.  4 
for  the  period  of  decaying  oscillations 
(oscillation  period  decreases  with  time),  and 
Fig.  5  for  RR  (RR  decreases  with  time), 

c)  As  seen  in  Fig.  2,  the  period  of  oscillations 
monotonically  increases  with  ^  for  the  second 
order  system.  On  the  other  hand,  as  seen  in 
Fig.  5  for  the  FOS,  starting  from  6.175  at  a  = 
1.34  the  duration  of  the  first  oscillations  first 
increases  and  it  reaches  to  its  maximum  value 
7.327  at  a  <  1.43  and  then  decreases  to  6.284 
at  a  =  1.99. 

V.  CONCLUSIONS 

Time  domain  step  response  characteristics  of  the  fractional 
order  system  with  unity  numerator  and  a  2-term 
denominator  polynomial  involving  a  single  fractional 
power  is  studied  in  this  paper.  Many  step  response 
characteristics  which  are  important  for  control 
engineering  such  as  rise  time,  settling  time,  delay  time, 
overshoot,  oscillation  period,  damping  time  constant  of 
oscillations,  and  some  others  are  investigated,  and  the 
results  are  presented  in  graphical  forms  by  figures.  The 
study  is  conducted  comparatively  by  considering  integer 
order  systems  of  1st  and  2nd  order  types.  It  is  shown  that 
the  same  simplicity  and  exp  licitness  present  for  second 
order  systems  do  not  exist  between  the  transfer  function 
parameters  and  the  step  response  characteristics  for  low  order 
fractional  order  control  systems.  The  results  bring  light  for 
the  design  such  systems  which  has  been  a  vacancy 
fulfilled  by  this  paper. 

The  derived  transcendental  equations  for  r  and  Tr  in  Eqs. 
(4a, 6a)  and  Eqs.  (4a, 6b),  respectively),  which  have,  as  far 
as  author  knowledge,  not  been  appeared  in  the  literature 
before,  are  derived  and  solved  numerically  to  obtain  the 
plots  for  relevant  characteristics  in  Fig.  2.  This  approach 
is  not  needed  for  obtaining  the  remaining  characteristics 
in  the  figure  for  the  necessary  formulas  are  presented 
explicitly.  For  the  simulation  of  the  FOS,  the 
subprograms  of  FOMCON  toolbox  of  Aleksei  Tepljakov 
[1]  integrated  with  MATLAB  R2017  [7]  is  used. 
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